On the approximation of Laplacian eigenvalues in graph disaggregation

Graph disaggregation is a technique used to address the high cost of computation for power law graphs on parallel processors. The few high-degree vertices are broken into multiple small-degree vertices, in order to allow for more efficient computation in parallel. In particular, we consider computations involving the graph Laplacian, which has significant applications, including diffusion mapping and graph partitioning, among others. We prove results regarding the spectral approximation of the Laplacian of the original graph by the Laplacian of the disaggregated graph. In addition, we construct an alternate disaggregation operator whose eigenvalues interlace those of the original Laplacian. Using this alternate operator, we construct a uniform preconditioner for the original graph Laplacian.     

Approximations of spectra of Schrödinger operators with complex potentials on ℝd

We study spectral approximations of Schrödinger operators T = ?Δ+Q with complex potentials on Ω = ?^d, or exterior domains Ω??^d, by domain truncation. Our weak assumptions cover wide classes of potentials Q for which T has discrete spectrum, of approximating domains Ω_n, and of boundary conditions on ?Ω_n such as mixed Dirichlet/Robin type. In particular, Re Q need not be bounded from below and Q may be singular. We prove generalized norm resolvent convergence and spectral exactness, i.e. approximation of all eigenvalues of T by those of the truncated operators T_n without spectral pollution. Moreover, we estimate the eigenvalue convergence rate and prove convergence of pseudospectra. Numerical computations for several examples, such as complex harmonic and cubic oscillators for d = 1,2,3, illustrate our results.     

An adaptive wavelet space‐time SUPG method for hyperbolic conservation laws

This article concerns with incorporating wavelet bases into existing streamline upwind Petrov‐Galerkin (SUPG) methods for the numerical solution of nonlinear hyperbolic conservation laws which are known to develop shock solutions. Here, we utilize an SUPG formulation using continuous Galerkin in space and discontinuous Galerkin in time. The main motivation for such a combination is that these methods have good stability properties thanks to adding diffusion in the direction of streamlines. But they are more expensive than explicit semidiscrete methods as they have to use space‐time formulations. Using wavelet bases we maintain the stability properties of SUPG methods while we reduce the cost of these methods significantly through natural adaptivity of wavelet expansions. In addition, wavelet bases have a hierarchical structure. We use this property to numerically investigate the hierarchical addition of an artificial diffusion for further stabilization in spirit of spectral diffusion. Furthermore, we add the hierarchical diffusion only in the vicinity of discontinuities using the feature of wavelet bases in detection of location of discontinuities. Also, we again use the last feature of the wavelet bases to perform a postprocessing using a denosing technique based on a minimization formulation to reduce Gibbs oscillations near discontinuities while keeping other regions intact. Finally, we show the performance of the proposed combination through some numerical examples including Burgers’, transport, and wave equations as well as systems of shallow water equations.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2062–2089, 2017     

Error and stability analysis of numerical solution for the time fractional nonlinear Schrödinger equation on scattered data of general‐shaped domains

In present work, a kind of spectral meshless radial point interpolation (SMRPI) technique is applied to the time fractional nonlinear Schrödinger equation in regular and irregular domains. The applied approach is based on erudite combination of meshless methods and spectral collocation techniques. The point interpolation method with the help of radial basis functions is used to construct shape functions which play as basis functions in the frame of SMRPI. It is proved the scheme is unconditionally stable with respect to the time variable in and also convergent by the order of convergence , . In the current work, the thin plate spline are used as the basis functions and to eliminate the nonlinearity, a simple predictor‐corrector (P‐C) scheme is performed. It is shown that the SMRPI solution, as a complex function, is suitable one for the time fractional nonlinear Schrödinger equation. The results of numerical experiments are compared to analytical solutions to confirm the reliable treatment of these stable solutions. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1043–1069, 2017     

Generalized Bessel functions: Theory and their applications

This paper presents 2 new classes of the Bessel functions on a compact domain [0,T] as generalized‐tempered Bessel functions of the first‐ and second‐kind which are denoted by GTBFs‐1 and GTBFs‐2. Two special cases corresponding to the GTBFs‐1 and GTBFs‐2 are considered. We first prove that these functions are as the solutions of 2 linear differential operators and then show that these operators are self‐adjoint on suitable domains. Some interesting properties of these sets of functions such as orthogonality, completeness, fractional derivatives and integrals, recursive relations, asymptotic formulas, and so on are proved in detail. Finally, these functions are performed to approximate some functions and also to solve 3 practical differential equations of fractionalorders.     

双调和算子特征值问题的混合三角谱元方法

本文针对双调和算子特征值问题设计了基于混合变分形式的三角谱元逼近格式,其基函数采用指标为(-1,-1,-1)的广义Koornwinder多项式.在H~1-及H_0~1-正交谱元投影的逼近理论基础上,我们建立了双调和算子特征值与特征函数的收敛性估计;它关于网格尺寸h是最优的,关于多项式次数M是次优的.然而,在H_0~2-正交谱元投影的最优估计假设前提下,关于M的次优收敛阶估计则提升为最优.此外,Koornwinder分片多项式逼近的结果还表明,在带权Besov空间范数的度量下,对于存在着区域角点奇性的双调和算子特征值问题,谱元方法的收敛阶能达到h-型有限元方法的2倍.最后,本文的数值实验结果展示了谱元逼近格式的高效性,同时也验证了相关理论的正确性.     

Asymptotic behavior of solutions to the compressible nematic liquid crystal system in ℝ3

In this paper, we study a nematic liquid crystals system in three-dimensional whole space ?³ and obtain the time decay rates of the higher-order spatial derivatives of the solution by the method of spectral analysis and energy estimates if the initial data belongs to L¹?³ additionally.     

一类弱奇性Volterra积分微分方程的级数展开数值解法

本文考察了一类弱奇性积分微分方程的级数展开数值解法,并给出了相应的收敛性分析.理论分析结果表明,若用已知函数的谱配置多项式逼近已知函数,那么方程的数值解以谱精度逼近方程的真解.数值实验数据也验证了这一理论分析结果.     

可修复人机系统的指数稳定性

研究了两相同部件温储备可修的人机系统,运用C_0半群的相关理论,对系统主算子的谱界进行估值.估算系统的算子产生的半群的增长界,然后运用了共尾的概念及相关的理论,得到了系统算子A+B的谱界与系统算子产生的半群的增长界相同.进而运用相关代数知识证得,0为系统算子的简单本征值,并分析了系统算子的谱分布,得到系统的指数稳定性.并研究了系统算子预解式的特性.对任意给定的δ0,γ=a+bi,-μ+δa_1≤a≤a_2,得到lim|b|→∞‖R(γ;A+B)‖=0.进而得到在~sRγ≥a_1的右半平面内相应于系统算子A+B的谱点由有限个本征值组成.